Does Bulk Viscosity Create a Viable Unified Dark Matter Model? 
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We investigate in detail the possibility that a single imperfect fluid with bulk viscosity can replace 
the need for separate dark matter and dark energy in cosmological models. With suitable choices 
of model parameters, we show that the background cosmology in this model can mimic that of a 
ACDM Universe to high precision. However, as the cosmic expansion goes through the decelerating- 
accelerating transition, the density perturbations in this fluid are rapidly damped out. We show 
that, although this does not significantly affect structure formation in baryonic matter, it makes the 
gravitational potential decay rapidly at late times, leading to modifications in predictions of cosmo- 
logical observables such as the CMB power spectrum and weak lensing. This model of unified dark 
matter is thus difficult to reconcile with astronomical observations. We also clarify the differences 
with respect to other unified dark matter models where the fluid is barotropic, i.e., p = p(p), such 
as the (generalized) Chaplygin gas model, and point out their observational difficulties. We also 
summarize the status of dark sector models with no new dynamical degrees of freedom introduced 
and discuss the problems with them. 

PACS numbers: 98.80-k, 95.36+x 



o 

i . 



> 
m 

(N 
O 
On 
O 



X 



I. INTRODUCTION 

In recent years the cosmological picture that over 95% 
of the energy in our universe is contributed by a dark 
sector has been supported independently by a number 
of observations, notably those of type la Supernovae 
(SN) luminosity distances, cosmic microwave background 
(CMB) anisotropics, the power spectrum of clustered 
matter, and weak lensing [H, 0j S 0]- This dark sec- 
tor can be further subdivided into dark matter and dark 
energy according to their different gravitational proper- 
ties. The concordance ACDM paradigm, in which dark 
matter is assumed to be weakly (or just gravitationally) 
interacting massive particles (the cold dark matter), and 
dark energy is a positive cosmological constant (A) or 
slowly varying scalar field, has been successful in con- 
fronting all these observational data sets. However, the 
smallness of the cosmological constant, and the fact that 
it only becomes dominant recently make this model con- 
ceptually unattractive and stimulates the examination of 
new models (for a review see Q). 

Since both cold dark matter (CDM) and dark energy 
are invisible and have as yet unknown origins, it is natural 
to consider the possibility that they are actually not two 
exotic matter species but just different aspects of a sin- 
gle fluid. These scenarios are frequently dubbed unified 
dark matter (UDM) models. In the context of general 
relativity they assume an equation of state p = p(p) or 
p = p(H), where H is the Hubble rate; in universes with 
zero spatial curvature these prescriptions are identical. A 
particular class of cosmology with this equation of state 
was investigated in the context of studies of cosmological 
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bulk viscosity [1, @, H, EH HH m which the viscosity 
coefficient is f](p) — ap m and the effective pressure in 
the Friedmann equation is p' = (7 — l)p — 3Hri(p), and 
of string production effects which mimic the effects 
of bulk viscosity of this form with m = 3/2 Q- 

The flat bulk viscous cosmologies also include as sub- 
cases the so called Chaplygin gas model {p = —Ap~ a ) 
, which is just a bulk viscosity for dust (7 = 1) 
with 771+1/2 = —a, and its generalizations to p+p = Bp x 
[l3|, which is just 7 = and m + 1/2 = A, or other 
functional forms p{p) [j| . The Chaplygin gas models are 
simple, with no new dynamical degrees of freedom, and 
yet produce an interesting time-evolution for the dark 
energy. However, it was shown subsequently that the 
density perturbation of this fluid will either blow up or 
experience rapid oscillatory damping at late times, so 
the models can be stringently constrained by the matter 
power spectrum [16[. This distinctive behaviour arises 
because there is a minimum possible total density of the 
universe at late times. 

There have been some explicit investigations in the 
bulk viscous cosmology in connection with the dark en- 
ergy problem. Fabris et al. [l8| investigated the case 
in which the viscosity coefficient has the form rj{p) = ap m 
and considered both background and perturbed evolu- 
tions of a universe dominated by viscous matter. How- 
ever, these studies were limited by certain simplifications. 
In Ref. for example, the authors obtained analyti- 
cal solution for p v (where the subscript v denotes viscous 
matter) as a function of scale factor a, under the assump- 
tion that no other matter species exist in the universe 
(the same assumption as [7|). In Ref. [l8| the authors 
added a baryonic matter species and estimated the cos- 
mological parameters in detail using Bayesian statistics; 
but their calculation was largely confined to the special 
case m = and used only supernovae data to derive the 
constraints. We believe that such an analysis needs to 
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be generalized. First, as will be seen below, m = is 
not necessarily the best fit for these models, even for the 
background cosmology. Second, no radiation component 
was included in previous analysis. This is a reasonable 
simplification as long as we are only concerned with the 
late-time cosmic expansion and only using the supernova 
data. But when observables are also related to the high- 
redshift features, such as the CMB shift parameter, the 
radiation must be taken into account. In fact, the prop- 
erties of the viscosity could significantly modify the early- 
time evolution of the universe if the late-time evolution is 
consistent with the supernova data, and so the CMB shift 
parameter could be effective in constraining the model. 
Third, a detailed study of the linear perturbations in the 
model is needed as a critical check of the model's fea- 
sibility. As we will see below, the UDM model based 
on viscous matter has several distinctive predictions re- 
garding the perturbations when compared to the ACDM 
paradigm or the Chaplygin gas model. Such a feature is 
very general, making the model barely compatible with 
observations. This indicates that a UDM model with- 
out new dynamical degrees of freedom is unlikely to be 
observationally acceptable. 

The rest of this paper is organized as follows. In ^ UTl 
we give the basic field equations which will be used in 
the subsequent analysis. In § IIIII we consider the back- 
ground evolution of a universe dominated by a viscous 
matter (in the presence of normal matter species such as 
baryons, photons and neutrinos). Using the SN data and 
CMB shift parameter, we find the best-fit model param- 
eters and show that these give a background cosmology 
very similar to the prediction of the ACDM paradigm. 
The § |IV] is devoted to the perturbed evolution of gen- 
eral bulk viscosity models, which we find to behave very 
differently from both ACDM and other UDM models. 
Finally, we discuss our results and conclude in § [V] Wc 
will frequently call the UDM fluid with p = p(p, H) 'the 
viscous dark matter'. 



II. THE FIELD EQUATIONS 

In this section we list the general field equations that 
govern the evolution of the cosmological background and 
its first-order perturbations in general relativity, which 
will be used in later sections. The perturbation equations 
will be given in the covariant and gauge invariant (CGI) 
formalism, using the method of 3 + 1 decomposition. 

The main idea of 3+ 1 decomposition is to make space- 
time splits of physical quantities with respect to the 4- 
velocity u a of an observer. The projection tensor h a b is 
defined as h a b = g ab — u a u b and can be used to obtain 
covariant tensors perpendicular to u. For example, the 
covariant spatial derivative V of a tensor field Tj,','.'" is 
defined as 

v" '/;,:::: ^ Kh) ■ ■ ■ hix, ■ ■ ■ h;vr-;t (1) 

The energy-momentum tensor and covariant derivative 



of the 4- velocity are decomposed respectively as 

Tab = 7r a6 + 2g( a u 6 ) + pu a Ub - ph ab , (2) 

V a U 6 = (T ab + KJab + -^0h ab + U a A b . (3) 

In the above, ir ab is the projected symmetric trace-free 
(PSTF) anisotropic stress, q a the vector heat flux vec- 
tor, p the isotropic pressure, a ab the PSTF shear ten- 
sor, zu ab — V[ a Uh], the vorticity, 9 = V c u c = 3a/ a (a is 
the mean expansion scale factor) the expansion scalar, 
and A b — ii b the acceleration; the overdot denotes time 
derivative expressed as <j) — u a \J a 4>, brackets mean an- 
tisymmetrisation, and parentheses symmetrization. The 
4- velocity normalization is chosen to be u a u a = 1- The 
quantities Tr a b, q a , p,p are referred to as dynamical quan- 
tities and a abl w ab ,9, A a as kinematical quantities. Note 
that the dynamical quantities can be obtained from the 
energy-momentum tensor T ab through the relations 

p = T ab U a u\ 

P = -\h ab T ab , 

q a = h d a u c T cd , 
ir ab = h c a hfT cd +ph ab . (4) 

Decomposing the Riemann tensor and making use the 
Einstein equations, we obtain, after linearization, a set 
of propagation and constraint equations governing the 
evolution of perturbed physical quantities. Here we shall 
only list the equations that will be used in later sections, 
and for more details we refer the reader to 

The first equation we will use is the Raychaudhuri 
equation 

0+\e 2 -VM + ~(p + 3p) = 0. (5) 

The second equation to be used involves the projected 
Ricci scalar R into the hypersurfaces orthogonal to u a , 
which can be expressed as 

2 o 

R = Inp- -6 2 . (6) 
3 

Since we are considering a spatially-flat universe, the spa- 
tial curvature vanishes (at background level) on large 
scales and so R = 0. Thus, from Eq. ©, we obtain 

\o 2 = k P) (7) 

which is the Friedmann equation that governs the expan- 
sion of the universe in standard general relativity. 

Furthermore, we will need the conservation equations 
for the energy-momentum tensor 

p+{p + p)9 + V a q a - 0, (8) 

qa+^Oq a + (p+p)A a -VaP + V b lTab = 0. (9) 
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Eqs. (0 [H [9]) involves both background and first-order 
perturbed quantities (such as A a , q a ). To obtain equa- 
tions for the background cosmology it is sufficient to ne- 
glect all first-order terms, while to obtain corresponding 
equations for the perturbation evolution we need to take 
the spatial derivatives V a [cf. Eq. ([1])] of these equations. 
For the perturbation analysis it is then more convenient 
to work in the k space because we shall confine ourselves 
to the linear regime where different fc-modes decouple. 
Following fl9[ ], we make the following harmonic expan- 
sions of our perturbation variables: 

c — ' a *■ — ' a 

k k 
k k 

VJ = J2^Q k a A a = Y^ k -AQ k a 

k k 

in which Q k is the eigenfunction of the comoving spatial 
Laplacian a 2 V 2 satisfying 

V 2 Q fc = %Q\ (10) 

or 

and Q k , Q k b are given by Q k = |V a Q k , Q k b = f V (o Q*. 
The perturbed version of Eq. ([8]) is 

Xa + {p + p){Z a - a6A a ) + (X a + X*)0 + aV Q V b g b =£1D) 

In § IIVI we shall use the harmonic expansion coefficients 
of Eqs. (0 [H [TTj) to derive the evolution equation for the 
density perturbation of the viscous dark matter. 

III. THE BACKGROUND EVOLUTION 

The field equations governing the background cosmic 
expansion are the Friedmann equation, the Raychaud- 
huri equation and the conservation equations for energy 
densities of the different matter species (baryons, radia- 
tion, and viscous dark matter). Not all these equations 
are independent, and we choose the Friedmann equation 
(here H — a/a — 8/5) 

3H 2 = 8ttG (pn + PB + Pr) , (12) 

and the conservation equations 

PB+3Hp B = 0, (13) 
p R + AHp R = 0, (14) 
pB + 3H(p D +p D ) = (15) 

as our starting point. In these equations pb,Pr and pd 
are respectively the energy densities of baryonic matter, 
radiation and the (viscous) dark matter. The pressure 
Pd = (7 — 1)/°d — 3aHp 7 £ is the effective pressure of the 
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FIG. 1: The joint constraints on m and /3 from Supernovae 
and CMB shift-parameter data. The dark grey and light grey 
areas denote the 68% and 95% confidence regions respectively. 
The current Hubble expansion rate Ho is marginalized over 
analytically and the other parameters used are J7r = 8.475 x 
10 -J and r^R + Q,b ~ 0.04 evaluated at the present time. 



dark matter, with 7, a, m being our free model parame- 
ters. 

We shall not follow [H[ by dividing pd into different 
components with different equation of states (EoS), al- 
though it can be useful mathematically, since that might 
hide the fact that there is only a single fluid with vary- 
ing EoS. When it comes to the perturbative evolution of 
this fluid, it will be misleading if one thinks that part of 
the fluid behaves as a cosmological constant which has no 
perturbation, while another part simply clusters as CDM, 
there being no interactions between them. We also note 
that the Hubble expansion rate H appears in the expres- 
sion of pd, which means that the EoS of the viscous dark 
matter depends also on the existence and properties of 
other matter species, and so we cannot neglect the effects 
of baryonic matter and radiation (at early times). 

Next, we estimate the free model parameters 7, a and 
m. We will set 7 = 1 in all the calculations below because 
we want the viscous dark matter to behave like CDM at 
(early) times when the correction term is not important. 
For to, earlier studies [?J showed that when to > 1/2 
the universe will start from a de Sitter phase and finally 
evolve towards power-law perfect-fluid dominated expan- 
sion, while if to < 1/2 it is just the opposite, with late- 
time approach to de Sitter evolution; the special to = 1/2 
case corresponds to power law evolution throughout. In 
fact, from the expression pn = —SaHp^ we can see 
that if pd dominates over other matter species (so that 

1 /2 

H cx p D ) then pn will be a constant if to = —1/2. We 
can estimate that the best fit value of to is around —1/2. 
In what follows, to will be taken as a free parameter to 
be constrained by data. Finally, a is obviously a dimen- 
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sional constant. To determine its value, we note that 
to explain the SNe data we need the candidate for dark 
energy to contribute an effective energy density and pres- 
sure of the same order as p\ in the ACDM model. Thus, 
we have 

Pdo — Pcdmo + PA, 
3aiJ PDo - PA 

in which the quantities on the left-hand sides are for our 
model while those on the right hand sides are for the 
ACDM model. A subscript o here denotes the present 
value of a quantity. Using the results Ocdm — 0.20 and 
51a — 0.76, we calculate from the above two equations 
that (1 = afloPoMo — 0-26. This gives us a sense about 
the magnitude of the dimcnsionlcss quantity (3, which is 
chosen as another model parameter instead of a and will 
be constrained below. 

With the above preliminaries, it is now straightforward 
to rewrite the conservation equation of the viscous dark 
matter, Eq. (fTS"]) , as 



Q + r bd e 



-3N 



-AN 



l + r b d + r rd 



1/2 



(16) 



where we have defined g = pdm/pdmo, Hd = Pbo/Pdmo, 
r r d = Pro/pdmo, and also used the Friedmann equation 
Eq. p2[) to substitute for the Hubble parameter, H . The 
star denotes a derivative with respect to N = log a. Be- 
cause there are in general no closed-form solutions to 
these cosmological equations, we will use Eq. (flB]) . to- 
gether with Eqs. (fl"2l [T51 [T4^) . in our subsequent numeri- 
cal calculation. Note that r b d and r rc i are constants which 
are fixed once Ob and Or are known (using the fact that 
the universe is spatially flat, so that Hdm = 1 — Ob — Or): 
^dm = 1 — Ob — Or): 



Or 



Or 



Tbd = 



1 - R - Or- 



Trd = 



1 - Or - T 



(17) 



A natural choice of the initial (final) condition of Eq. JTBJ 
is g(N = 0) = 1. 

We have used the supernovae luminosity distance data 
Q and CMB shift parameter [20] to constrain the model 
parameters m and (3, analytically marginalizing the cur- 
rent Hubble expansion rate Hq and assuming the baryon 
density today as fixed by BBN and measurements of 
the light element abundances. The result is shown in 
Fig. [T] and the best-fit parameters we find are (m, (3) = 
(—0.4, 0.236) which lie close to our estimate above (m = 
—0.5, j3 = 0.264). Fig. [5] shows the cosmic evolution of 
the fractional energy densities in the bulk viscous model 
with the above best-fit parameters It can be seen there 
that the viscous model mimics the concordance ACDM 
paradigm extremely well all through its cosmic history. 
The model therefore appears to work well as a description 
of the background cosmology. 




FIG. 2: (Color Online) Solid lines: The evolution of the frac- 
tional energy densities of viscous dark matter Qd (black), 
baryons Qb (red), and radiation (including photons and mass- 
less neutrinos) Qr (blue), versus 1 + z where z is the redshift. 
The model parameters here are chosen as m = —0.4 and 
f3 = 0.236. Dashed lines: the same evolutions for the concor- 
dance ACDM model. Here, fin denotes the fractional energy 
density of the dark sector, i.e., dark energy (a cosmological 
constant) plus cold dark matter. The other parameters used 
(for both models) are Q R = 8.475 x 10" 5 and fi a + fi B = 0.04, 
evaluated at the present time. 



IV. THE EVOLUTION OF FIRST-ORDER 
DENSITY PERTURBATIONS 

Despite the excellent coincidence between the viscous 
dark matter and ACDM models for the background evo- 
lution found in the last section, we should also investigate 
the formation of large-scale structure to test whether the 
viscous model is also a feasible model for dark energy. In 
this section we show that, generally, a bulk viscosity de- 
pending on the energy density, i.e., pu = — ??(pD)V a u a , 
will significantly influence the formation and evolution 
of large-scale cosmological structure. This very differ- 
ent prediction from that of the ACDM model indicates 
that stringent constraints can be placed on the viability 
of the bulk viscosity models using observational data on 
the CMB spectrum, matter power spectrum, and weak 
gravitational lensing. 

Let us first concentrate on the special case consid- 
ered above, with pu = —ap^\7 a u a — —3aHp^. For 
an observer comoving with dark matter particles (with 
4- velocity u a ), the energy momentum tensor could be 
written as 



Tab = PDU a U b - PD(g a b - U a U b ). 



(18) 



In the ACDM paradigm, if one chooses the observer to be 
comoving with the dark matter particles as above, then 
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FIG. 3: (Color Online) Upper panel: The evolution of the den- 
sity perturbation Ad for the viscous dark matter (the solid 
curves) , as compared to the predictions for the cold dark mat- 
ter model (the dashed curves). The results for 4 different 
length scales (k = 0.005,0.01,0.05,0.1 Mpc" 1 ) are displayed 
as labelled beside the curves. The model parameters are cho- 
sen to be (m = —0.4, (3 = 0.236). Lower panel: The same 
but for the baryon density perturbation. 



obviously the peculiar velocity is zero, vd — 0, which im- 
plies, by the conservation of energy- momentum tensor, 
that A a = for this observer, and there is no accelera- 
tion. In the bulk viscous model, however, vd — and 
A a = are different choices of frame which can be used 
in numerical calculations. Here, for convenience, we will 
choose the vd = frame, in which, again from the con- 
servation of energy, A a ^ 0, which is easy to understand 
because the dark-matter particles themselves have inter- 
actions (the pressure pd) and cannot be acceleration- free. 

Taking the spatial derivative of pu and picking the 
harmonic coefficients, we obtain 



XI = 



mp D A D - ap^-Z 
a 



where H = a' /a with ' = djdr (r is the conformal time 
defined by adr = dt) and Ap = X^/p^ is the density 
contrast of the viscous dark matter. 

As the anisotropic stress vanishes up to first order in 
perturbations, from Eq. ([9]) we obtain 



(p D + PD )A = XP. 
Similarly, from Eq. (fTTj) we have 

Ap + (1 + w)kZ - 3wHA L 







(20) 



(21) 



where we have defined the zero-order EoS parameter w 
w(a) = pd/pd for the viscous dark matter so that 



-3a Hp^ 



(22) 



Finally, taking the spatial derivative of the Raychaudhuri 
equation Eq. ((5]), we get 



kZ' + KHZ - k 2 A + 3(H' - H 2 )A 



'-{X + 3X p )al23) 



We shall assume here that the universe is dominated 
by the viscous dark matter, so that on the right-hand side 
of Eq. (f2"5)) X and X p can be replaced by X^> and X^, 
respectively. In more general cases it is straightforward 
to include contributions from other matter species. Then, 
from Eqs. 1(15]). (j2"D|) , (|2T|) and (J2S]), we can eliminate Z 
and A to obtain 

Ad + [C'i(a) + k 2 C 2 (a)} A D + [C 3 (a) + k 2 Ci(a)] A D ^24) 
where we have defined 

-,2 



Ci(o) 

C 2 (a) 
C 3 (a) 



H 



npua" 
2H 



n 



H 



f 



- w 
w 



— SwTL; 
1 



3(1 + 1/;) W 



(25) 
(26) 



npiya 



[(1 + w)(l + 3mw) + 3w 2 
—H + wT-A - 3wH 2 



-3w(H' -H 2 )[m 



w 



Ci(a) — w I m + 



1 



w 



(27) 
(28) 



It is straightforward to check that when w = (whether 
or not m = 0), this reduces to the evolution equation 
for the CDM density contrast, which does not depend 
on k, indicating that the evolution of the CDM den- 
sity contrast is scale-independent. For the viscous dark 
matter, however, we see that the equation depends on k 
in the coefficients of both A' and A. Eq. is effec- 
tively an equation for a damped simple oscillator, with 
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time-varying and scale-dependent frequency and damp- 
ing force, and no driving force. On small scales (k 3> TL) 
and at late times (when w ~ 0{— 1)), the C2(a) and 
04(a) terms dominate the coefficients of A' and A, re- 
spectively; it is then easy to show that the oscillator is 
over-damped so that its amplitude decays to zero rapidly 
with no oscillations. This qualitative feature can be seen 
in the upper panel of Fig. [31 where we plot the evolution 
of the dark matter density contrast, Ad, for four differ- 
ent length scales, k = 0.005, 0.01, 0.05, 0.1 Mpc -1 , for our 
best- fit model parameters (m = —0.4, (3 — 0.236). We see 
that the evolution of Ad deviates from that predicted 
by the ACDM model only at late times, but nonetheless 
significantly. The scale-dependent evolution of density 
contrast is actually a general feature in the UDM models 
without new dynamical degrees of freedom [l6|, [U |22| 

The lower panel of Fig. [3] displays the evolution of the 
baryon density contrast Ab- Because viscous dark mat- 
ter clusters more weakly than cold dark matter, the grav- 
itational potential that the baryons lie in is also weaker, 
making the baryons less clustered. However, since this 
effect is indirect, the deviation of Ab from the ACDM 
prediction is considerably smaller than that of Ad- Note 
that the galaxy surveys actually measure the clustering 
of luminous (baryonic) matter, rather than that of dark 
matter. Consequently, these measurements can only be 
applied to Ab, which is just weakly dependent on the 
model parameters (23[. Nevertheless, as in the bulk vis- 
cous model, both the amplitude and the shape of bary- 
onic matter power spectrum are distinct from those pre- 
dicted by the ACDM paradigm, and therefore we expect 
that stringent constraints can be obtained from it. 

The analysis above is generalizable to an arbitrary 
choice of the viscosity function i](pv). To do this, we 
just need to define 



rh(a) 



Pu dr)(p u ) 
v(pd) dp D 



(29) 



and replace the constants m in Eqs. (f2"Tl |2"5)) with rh(a). 
Then, Eq. (|24|) describes the evolution of the dark matter 
density contrast in the general case of 77(^0)- 

There are a couple of points to be noted about Eq. (f2"4")l . 
First, it is clear that m = does not guarantee that Ci 
and C4 are zero, and the EoS of the viscous dark matter 
w is important. Since w =/= 0, if the viscous dark matter 
is responsible for accelerating the expansion of the uni- 
verse, our conclusion that the evolution of dark matter 
density contrast will be sensitively scale-dependent and 
deviate significantly from ACDM will hold in general. 
Note that this is different from the analysis of Ref. [18l |. 
which considers the special case of m = and concludes 
that the matter power spectrum in this model behaves 
well. Second, Eq. (f2~4"|) is qualitatively different from its 
counterparts in other models which unify dark matter 
and dark energy, such as the (generalized) Chaplygin gas 
model. In the latter we have generally C2 = 0, and so 
the equation describes an under-damped oscillator rather 
than an over-damped one. Consequently, the dark matter 
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FIG. 4: (Color Online) The evolution of the gravitational po- 
tential <f> in the bulk viscosity model (solid curve) as compared 
to the results for the standard ACDM paradigm (dashed 
curves). The results for 4 different length scales (k = 0.0005, 
0.001,0.005 and 0.01 Mpc" 1 ) are shown. Clearly <j> decays 
much faster in the bulk viscosity model. 
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FIG. 5: The CMB temperature fluctuation spectrum of the 
bulk viscosity model (the solid curve) as compared to the 
predictions of the ACDM paradigm (the dashed curve). The 
model parameters are chosen to be m = —0.4 and /3 = 0.236, 
and all other parameters are the same in the two models. 



density contrast oscillates with (rapidly) decaying ampli- 
tude (or blows up if C4 is negative), in contrast to the 
monotonic decay in our model here. Tracing the deriva- 
tion of Eq. (|2"4"]l , it is easy to find that the C2 term comes 
from the term with Z in Eq. (|19[) . which is created by 
the fact that pd oc V a u a . We also note that Eqs. 
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[26)) are independent of m and hence independent of the 
functional form of rj(pr,). As long as C4 > 0, we will 
obtain a qualitative picture similar to the one given in 
Fig.H 

Although the very different evolution of Ad from that 
in ACDM is not directly reflected in the observed galaxy 
power spectrum, it will modify the gravitational potential 
and subsequently affect observables such as the CMB, 
weak lensing and CMB-galaxy cross-correlation. As the 
viscous dark matter contributes the majority of the total 
energy in the universe, a decay in its density contrast as 
in Fig. [3] will also drive the gravitational potential <f> to 
decay significantly. This point is verified in Fig. 01 which 
clearly shows that the decay of 4> is much faster than that 
in the ACDM model. 

The fast decay of <f> will enhance the integrated Sachs- 
Wolfe (1SW) effect, which then contributes a source term 
cx J T ° <f>' 'je[k(To — r)]dr to the CMB fluctuations, where 
ji(kr) is the spherical Bessel function and </>' the (con- 
formal) time derivative of <f>. This corresponds to more 
power in the CMB spectrum on large scales [24j . as dis- 
played in Fig. [5l Note that the positions of the acoustic 
peaks for the two models are the same in Fig. [5J be- 
cause their background evolutions are almost indistin- 
guishable. We can see that comparison with CMB data 
could provide strong restrictions on the present model 
as well. Again, because the damping in Ad is a general 
feature in UDM models based on bulk viscosity, so are 
the enhancements of the ISW effect and the low-£ CMB 
power. In principle, the weak-lensing convergence power 
spectrum, which reflects the (projected) potential distri- 
bution along the line of sight, will also be modified (as 
compared to the ACDM result) significantly by the rapid 
decay of the potential. This is not considered here be- 
cause the CMB spectrum itself already places stringent 
constraint on the model. 



V. DISCUSSION AND CONCLUSIONS 

To summarize, in this work we have investigated the 
background cosmological evolution and large-scale struc- 
ture formation in the bulk viscous models designed as 
alternatives to dark energy. A bulk viscosity generates 
an effective pressure p = —3rj(p)H which, when the func- 
tion rj{p) is appropriately chosen, could drive the acceler- 
ation of the cosmic expansion. Our numerical calculation 
focuses on a particular choice rj(p) = ap m with en, m be- 
ing constants. Such a choice has been considered before 
in the literature Q , in the contexts of both inflationary 
and late-time accelerating universes. It is interesting to 
note that with suitable values of m, the viscous mat- 
ter behaves as cold dark matter at early times and as a 
cosmological constant in the future; thus this model nat- 
urally unifies dark matter and dark energy, at least at 
the background level. 

We have used the measurements of supernovae lumi- 
nosity distance and CMB shift parameter to constrain 



the model parameters and found that for the best-fitting 
parameters, m = —0.4 and (3 — 0.236 (/3 <x a), the back- 
ground evolution of the universe is almost the same as 
that predicted by the ACDM model. From this view- 
point it seems that the model is a feasible alternative to 
an explicit cosmological constant. 

However, when it comes to the linear perturbations, 
the viscous dark matter starts to behave very differently 
from cold dark matter. The pressure of the viscous dark 
matter tends to resist the growth of the density contrast, 
and when it becomes significant (at late times) this effect 
can rapidly damp out the density perturbations, particu- 
larly on small scales. Even though this smoothing of the 
viscous dark matter density perturbation cannot be seen 
directly, it could reduce the growth rate of density pertur- 
bations in the luminous matter (which can be measured 
by the galaxy power spectrum) and drive the fast decay 
of the gravitational potential fluctuations, which subse- 
quently modifies the large-scale CMB spectrum, weak 
lensing and CMB-galaxy cross-correlations. 

Note that the model we consider has no explicit ACDM 
limit, i.e., one cannot adjust the model parameters a and 
m to make the model reduce to ACDM exactly, unless a 
cosmological constant is added and the limit a — > is 
taken. The latter case, in which we have both an explicit 
cosmological constant and viscous matter, is not particu- 
larly appealing because it will introduce more complexity 
without solving any of the problems of ACDM. Rather, 
we are interested in whether the viscous dark matter 
alone could replace ACDM completely. This means that 
the EoS parameter w will necessarily be of order —1 at 
late times. According to our analysis in § HVI [cf. Eqs. [|24l 
[2"8])]. it is the value of w that determines the evolution 
of Ad, and so we expect that all the qualitative pictures 
given in § IIVI will remain in place in any attempts to 
replace dark matter and dark energy with bulk viscous 
matter alone. Viscous matter is therefore not a successful 
contender for the dark sector. 

The bulk viscosity model is another candidate to ex- 
plain the dark energy without introducing dynamical de- 
grees of freedom. In order to explain the accelerating cos- 
mic expansion, one needs a negative (effective) pressure. 
If no new dynamical degree of freedom is introduced, the 
negative pressure should be either a constant, or a func- 
tion of /Ode (the energy density of the newly added mat- 
ter), p m (the energy density of the existing matter) or 
H, or combinations of these variables, which are all the 
possible variables in the background cosmology [25| (see 
[26l |27| for a counterexample however, where a new de- 
gree of freedom is added but is made non-dynamical). 
Hence, we have the following conclusions for each of the 
allowed prescriptions for a form of newly added matter 

1. The case p = const, is simply a cosmological con- 
stant. 

2. An example of the case p = jp(ppE) is the (gener- 
alized) Chaplygin gas model [14|, UM- The large- 
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scale structure formation here has been considered 
in [l6|, [U, [22|]. In this class of models we have 
X p oc X in general, and this leads to a term propor- 
tional to k 2 A in the evolution equation for the den- 
sity contrast A of the newly added matter, and this 
term dominates on small scales (k S> 7i). Depend- 
ing on the sign of this new term, A will either blow 
up or oscillate and decay rapidly on small scales. 

Note that this class of models are frequently consid- 
ered as either UDM models or dark energy models. 
In the former case, as shown in (l6| . the constraints 
on the models are particularly stringent (note how- 
ever that, as pointed out in 31, 32], in these models 
the nonlinear corrections can be important, which 
will make the simple linear treatments inaccurate 
or even incorrect, and will potentially significantly 
modify the background evolution as well, depend- 
ing on the parameters used). In the latter case, den- 
sity perturbations of CDM and baryons may not be 
affected significantly, much like the fact that in the 
above viscous model the baryons are not greatly 
affected. However, there generally will be other 
distinct new features of the model [21]. 

3. Examples for the case p = p(p m ) include the Car- 
dassian model 33], the Palatini f(R) gravity [34j ] 
and the u> = —3/2 Brans-Dicke theory with a po- 
tential, the linear perturbations of which have been 
considered in [H, [36|, H3, [M, UK . In these models, 
by a similar argument as in the above, there will 
be a term proportional to k 2 A m in the evolution 
equation for the existing matter density perturba- 
tion A m . As a result, A m will experience blowing- 
up or rapidly decaying oscillation on small scales at 
late times. Note, however, that in this class of mod- 
els the averaging over the microscopic structure of 
the (existing) matter distribution may be a serious 
issue, rendering the appropriately averaged cosmo- 
logical behaviour very different from naive predic- 
tions [4(| Hl| (in a way similar t o |32| though with 
some differences, namely in [|(| |41| the averaging 
is over microscopic scales while in [32| it is over as- 
tronomical scales) . This is because p for the newly 
added matter depends algebraically on the density 
of normal matter particles, and could be very differ- 
ent inside and outside the distribution of the latter. 



As the normal matter particles only occupy a tiny 
portion of the total volume of the universe, after av- 
eraging, p should be dominated by its value outside 
the normal matter particles (i.e., in the vacuum), 
which is likely to be a constant. 

4. The case p — p(H) includes bulk viscosit y m odel 
with 7](p) = const, considered in [l8|, l42l. |43|. |44| . 
which is a special case for the general model we 
considered in this work, with p = p(pBE, H). Here, 
the dependence of p on H results in the evolution 
equation of A for the newly added matter acquiring 
both a term proportional to k 2 A and one propor- 
tional to k 2 A'. Consequently, on small scales at 
late times, A will decay rapidly without oscillation 
like an over-damped oscillator. We have shown that 
the UDM model based on this is generally unable to 
pass several cosmological tests. Note that as in the 
case of p = p(pde), the averaging issue again needs 
to be taken into account if more precise predictions 
are to be obtained [3l|, [HJ , though its full effect and 
significance need careful nonlinear studies such as 
N-body simulations. 

In all such cases (except p — const.), we have seen 
that the evolution of A (for either the existing matter or 
the newly added matter) becomes very irregular. If the 
matter with irregular perturbation evolution makes a sig- 
nificant contribution to the total energy of the universe, 
then this will lead to large modifications to the ACDM 
predictions for the matter power spectrum, CMB, weak 
lensing and similar tests. This indicates that a viable dy- 
namical UDM model is likely to involve extra dynamical 
degrees of freedom in contrast to that provided by a (gen- 
eralized) Chaplygin gas or bulk viscosity |45l |46|. l47l |48| . 
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